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Spectroscopic and photometric evidence indicates that Type la supernovae (SNe la) are the thermonu¬ 
clear explosions of accreting white dwarfs. However, the progenitor binary systems and hydrodynamical 
models for SNe la are still controversial. The relatively uniform light curves and spectral evolution of SNe 
la have led to their use as a standard candle for determining cosmological parameters, such as the Hubble 
constant, the density parameter, and the cosmological constant. Recent progress includes the calibration 
of the absolute maximum brightness of SNe la with the Hubble Space Telescope, the reduction of the 
dispersion in the Hubble diagram through the use of the relation between the light curve shape and the 
maximum brightness of SNe la, and the discovery of many SNe la with high red shifts. 


Supernovae are classified spectroscopically as Type I if 
they have no hydrogen lines in their optical spectra and 
Type H if they have hydrogen lines in their optical spec¬ 
tra. Type I supernovae (SNe I) are further subclassified 
into types la, Ib, and Ic on the basis of spectra observed 
early in their explosion (early-time spectra) (|^). SNe 
la are characterized by the presence of a deep Si H ab¬ 
sorption line near wavelength 6150 A(Fig. 1), and their 
late-time spectra are dominated by strong blends of Fe 
emission lines. SNe Ib and Ic, in contrast, do not show 
this Si line. Moderately strong He I lines, especially at 
5876 A, distinguish SNe Ib from SNe Ic in early-time 
spectra; that is, SNe Ib exhibit absorption lines of He 
I, whereas these lines are weak or absent in SNe Ic (P. 
SNe H, Ib, and Ic are now generally thought to result 
from the explosion of massive stars—SNe H from single 
stars and SNe Ib and Ic from binary stars. There are 
spectroscopic and photometric indications that SNe la 
originate from white dwarfs that are composed of C -I- 
O with strongly degenerate electrons and have accreted 
sufficient mass from a companion to trigger an explosion. 

However, the progenitor systems and hydrodynami¬ 
cal models for SNe la are still controversial. Many is¬ 
sues need to be resolved including (i) double degener- 
ate(DD) versus single degenerate(SD) scenarios, that is, 
whether the companion of the white dwarf is also an 
electron-degenerate white dwarf or is a non-degenerate 
[main sequence (MS) or evolved red giant] star; (ii) Chan¬ 
drasekhar mass (Ch) (H) versus sub-Chandrasekhar mass 
(sub-Ch) models; and (iii) the explosion mechanism in 
the Ch models. The answers to these questions could 
lead to improved understanding of the complicated evo¬ 
lution of close binaries as well as the physics of thermonu¬ 
clear explosion (P. The relatively uniform light curves 
and spectral evolution of SNe la have led to their use 


as a standard candle to determine the Hubble constant 
(iJo), the density parameter (Hm), and the cosmological 
constant (Ha) (|J). SNe la are important standard can¬ 
dles because they are bright enough to be observed out 
to a redshift (z) of ^ 1. Variations of light curves and 
spectra among SNe la have recently received attention. 
SNe la with higher maximum brightnesses tend to show 
a slower decline in their light curves (||, ||; Fig. 2). This 
review summarizes our current understanding of SNe la 
and focuses on these controversial issues. 

Progenitors 

White dwarfs composed of C-l-0 are formed from inter¬ 
mediate mass stars (M < 8 M©, where Mq is the mass 
of our sun), undergo cooling, and eventually become dark 
matter, as they evolve towards fainter luminosities. In a 
close binary system the white dwarf evolves differently 
because the companion star expands to transfer matter 
to the white dwarf; the accreting white dwarfs are re¬ 
juvenated and in certain cases undergo thermonuclear 
explosions to give rise to SNe la. Theoretically the Ch 
white dwarf models and the sub-Ch models have been 
considered to explain the origin of SNe la (|^). Various 
evolutionary scenarios have been proposed, including (i) 
merging of double C-l-0 white dwarfs with a combined 
mass exceeding the Ch limit (a DD scenario) (H) and 
(ii) accretion of H or He by mass transfer from a binary 
companion at a relatively high rate (an SD scenario) (|^, 
!)■ 

The DD-Ch scenario is favored from the theoretical es¬ 
timate of frequencies of the occurrence of SNe la. This 
has stimulated a search for DD systems, but few such 
systems have been discovered, and their combined mass 
is smaller than Ch (^. For an SD scenario, the accret- 
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ing white dwarf undergoes H burning near the surface, 
which increases or decreases the white dwarf mass de¬ 
pending mainly on the accretion rate, M (^. Ch mass 
white dwarfs can be formed with a relatively high ac¬ 
cretion rate such as M w 10“® to 10“® Mq year“^ 
because of relatively small mass ejection after H or He 
shell burning. Recently, it was found that this range of 
M can be extended to a rate faster than 10“® Mq year“^ 
With such a high rate, the white dwarf undergoes 
wind mass loss without expanding its radius, thereby in¬ 
creasing its mass by steady H burning. The SD-sub-Ch 
scenario—that is, the explosion of sub-Ch white dwarfs— 
is a possible outcome of accretion with 4 x 10“® Mq 
yr“^ ^M Mq year"In this case, the ignited 

He shell flash is strong enough to initiate an off-center 
He detonation (0 12), which induces various types of 
explosions. 

Promising candidates for the white dwarf’s companion 
stars for the SD scenario are (i) stars with MS masses of 
1 to 1.5 Mq, which fill the Roche lobe when they evolve 
through red-giants, and (ii) the 2 to 3 Mq stars, which 
undergo mass transfer near the MS. These cases of rela¬ 
tively fast accretion onto the white dwarf can correspond 
to super-soft x-ray sources ( 0 ), symbiotic stars ( 0 ): or 
both. 

To discriminate among the SD-sub-Ch, SD-Ch, and 
DD-Ch scenarios, we used photometric and spectroscopic 
diagnostics. If any H or He is detected, the DD model 
could be ruled out. Although high-velocity H has not 
been observed in any SNe la, the upper limit to the H 
abundance ( 10 ^ Mq) is still too high to rule out 

the SD models. The H-rich materials in the companion 
star may be engulfed in the exploding material; thus the 
detection of low-velocity H could be critical 

Explosion Models 


Once explosive nuclear burning is ignited, it induces ther¬ 
monuclear explosion of the white dwarf. The outcome of 
this explosion depends on how the nuclear flames and/or 
shock waves propagate for the Ch and sub-Ch mass white 
dwarfs. The physics involved in these processes is rather 
complex, and multi-dimensional simulations have been 
conducted to understand these processes The re¬ 

sults are still preliminary, and so we focus here on those 
models that can account for the basic features of SNe la, 
namely, an explosion energy of ~ 10^'^ J, the synthe¬ 
sis of a large amount of ®®Ni (0.4 to 1 Mq), and the 
production of a substantial amount of intermediate-mass 
elements at expansion velocities of ~ 10,000 km near 
the maximum brightness of the SN explosion. 
Chandrasekhar mass models: Carbon is ignited in the 
central region of the white dwarf when the central den¬ 
sity exceeds ~ 10® g cm"®. Because of strong electron 


degeneracy, C burning is so explosive that it incinerates 
the material into Fe peak elements Afterwards, the 
explosive nuclear flame propagates outwards. The flame 
front is subject to various types of instabilities, including 
thermal instabilities, the Landau-Darrius instability, the 
Rayleigh-Taylor instability, and the Kelvin-Helmholtz in¬ 
stability (0)- The flame speed depends on the develop¬ 
ment of these instabilities and the resulting turbulence, 
requiring further extensive simulations on large and small 
scales as well as suitable modeling of the turbulence. Be¬ 
hind the flame front, materials undergo explosive nuclear 
burning of Si, O, Ne, and C. The nucleosynthesis prod¬ 
ucts depend mainly on peak temperatures, which in turn 
depend on the densities encountered by the flame. For 
densities from 10^® to 10® g cm"®, the products range 
from Fe-peak elements (mostly ®®Ni) to intermediate- 
mass elements (Ca, Ar, S, and Si), 0-Ne-Mg, and C-l-0. 

During the subsonic propagation of the deflagration 
wave, the densities of the whole white dwarf are de¬ 
creasing because of expansion. Therefore, the densities 
encountered by the flame are determined by its speed, 
which is still uncertain. Several plausible models have 
been presented with one-dimensional codes (0. 12). In 
carbon deflagration models, such as model W7 (0), the 
average flame speed is as high as one-fifth of the sound 
speed. A sequence of nucleosynthesis reactions produces 
®®Ni, Ca-S-Si, 0-Ne-Mg, and C-l-0 behind the deflagra¬ 
tion wave. Because of fast propagation of the flame, the 
transition to the detonation could be induced by shock 
compression in the outer low density layers, which would 
produce some variations of ®®Ni mass and distribution, 
as seen in peculiar SNe la such as SN1991T (Fig.l). 
In contrast, delayed detonation models assume that the 
early propagation of the deflagration is as slow as a few 
percent of the sound speed, producing ~ 0.1 Mq of ®®Ni, 
and hence the transition from deflagration to detonation 
could occur at densities of 1 x 10^ to 3 x 10^ g cm"®. 
In this case, the bulk of the white dwarf has expanded to 
lower densities so that the detonation wave synthesizes 
Fe-peak elements and intermediate mass elements &■ 
In the pulsating delayed detonation model, the transition 
to detonation is assumed to occur near the maximum 
compression (|^). Possible variations of the transition 
density could produce variations of the ®®Ni mass pro¬ 
duced. 

Sub-Chandrasekhar mass models: A central C detonation 
can be initiated by a He detonation-induced shock wave 
(12, 21). Merging of double white dwarfs could ignite a 
central C detonation ( p^ . If the white dwarf mass is less 
than 1.07 Mq, the densities of the white dwarf matter 
encountered by the detonation wave can be suitable to 
produce sufficient amounts of ®®Ni and Si-Ca-0. Varia¬ 
tions of the white dwarf mass could cause the variations 
of the light curves of the resultant SN la. 

Constraints on still uncertain parameters in these mod- 
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els, such as the central ignition density, the flame speed, 
and the deflagration-detonation transition density, can 
be provided by comparisons of theoretical spectra and 
light curves with observations and by comparisons of nu¬ 
cleosynthesis consequences with solar isotopic ratios. A 
combination of processes from various models may be 
responsible for the variations in characteristics of the ob¬ 
served SNe la. 

Spectra and Light Curves 

Spectra: Spectra of SNe la are generally homogeneous 
but show some important variations. Because SNe la do 
not have a thick H-rich envelope, elements newly synthe¬ 
sized during the explosion can be observed in the spectra. 
Thus, comparison between the synthetic spectra and ob¬ 
servations is a powerful diagnostic of the dynamics and 
nucleosynthesis suggested by the models. Non-local ther¬ 
modynamical equilibrium(LTE) spectra have been calcu¬ 
lated and can be used to compare with observations. A 
typical example (Fig. 3) shows agreement between the C 
deflagration model (unmixed W7) and the observed opti¬ 
cal spectra of SNe 1992A and 1994D (^^. The material 
velocity at the stellar surface near maximum light is ~ 
10,000 km s“^and the spectral features are identified as 
those of Fe, Ca, S, Si, Mg, and O. This implies that the 
abundance distribution in the velocity space may be sim¬ 
ilar to W7. The synthetic spectra for the sub-Ch models 
seem to be less satisfactory ( p^ pj). 

Regarding the heterogeneity in the spectra, SNe 1991T 
(Fig. 1) and 1991bg are the two extreme examples that 
have clearly revealed the presence of spectroscopically 
peculiar SNe la (|^). The pre-maximum light spectra 
of SNe la show a significant variation of the composition 
and expansion velocities of the outermost layers, whereas 
the post-maximum spectra are relatively uniform except 
for SN 1991bg. The different ®®Ni mass may produce 
the variation of the spectra (p^. Further analysis could 
provide the abundance distribution of various elements 
(Fe, Ca, Si, Mg, O, C, and others) in velocity space (p^. 

Light curves: Recent high-quality charge-coupled device 
(CCD) observations of SNe la have, on one hand, sup¬ 
ported the basic homogeneity of the optical light curve 
shape of SNe la. Good examples include SNe 1980N 
and 1981D in the galaxy NGC 1316, whose light curve 
shapes and brightnesses were almost identical. On the 
other hand, the observations have established significant 
variations of the maximum brightness, the light curve 
shape, and their correlation (Fig.2). In theoretical mod¬ 
els, the explosion energy goes into the kinetic energy of 
expansion, E. The light curves are powered by the ra¬ 
dioactive decay sequence ^®Ni —> ®®Co —> ^®Fe. The cal¬ 
culated light curve reaches its peak at about 15 to 20 
days after the explosion and declines because of the in¬ 


creasing transparency of the ejecta to gamma rays as 
well as the decreasing input of radioactivity. The light 
curve shape depends mainly on the effective diffusion 
time Tm oc (tt;M/uexpc)^^^, where k is opacity, c is the 
speed of light, and Uexp oc (|l^. For longer Tm, 

the decline of the light curve is slower. For larger Mni, 
the brightness is higher. For the Ch models, typical val¬ 
ues of Mni and E are similar to those of W7, that is, Mni 
= 0.6 ± 0.1 Mq and F = 1.3 ± 0.1 x lO^^ J. 

The observed uniformity in the light curve implies a 
uniformity of M, A, and Mni. A homogeneous M and 
E can be naturally accounted for with the Ch models, 
although the central density can vary with a small differ¬ 
ence of mass near the Ch limit. For the sub-Ch models, 
the accretion process must somehow choose a relatively 
narrow mass range, such as M « 1.0 to 1.1 Mq. The 
reported variation of the maximum brightness of SNe la 
®|) may be a result of the variations of the ^®Ni mass. 
The amount of ^®Ni depends on (i) the flame speed in 
the fast-deflagration models, (ii) the location and den¬ 
sity of the transition to a detonation in the delayed or 
late detonation models, and (iii) the white dwarf mass 
in the sub-Ch models. Because of variations of these pa¬ 
rameters, a variation of the ®®Ni mass may be possible. 
For the Ch models, the dependence of optical opacities 
on the temperature is important (|^) . A smaller mass of 
®®Ni means less heating of the surface layer, thus lowering 
K and Tm ■ This effect might explain the relation between 
brightness and the decline rate of the light curve. For 
the sub-Ch models, Mni varies approximately in propor¬ 
tion to M. Then for larger M, Mni is larger and Tm is 
longer; in other words, brighter SNe la tend to decline 
more slowly (p^. 

Nucleosynthesis in Galactic 
Chemical Evolution 

Supernovae of different types have different progenitors, 
thus producing different heavy elements on different time 
scales during the chemical evolution of galaxies. A rea¬ 
sonable mixture of the heavy-element yields from SNe 
la and SNe II should be able to explain the solar abun¬ 
dance pattern of heavy elements from O to the Fe group. 
Nucleosynthesis products of SNe la and SNe II can be 
combined with various ratios and compared with solar 
abundances of heavy elements and their isotopes. If the 
nucleosynthesis products of SNe II as a function of stellar 
masses ( ^^ are adopted with an upper mass limit of 50 
M 0 and the nucleosynthesis products of SNe la are those 
from model W7, the best fit to known solar abundances 
is obtained if the number of SNe la that have occurred 
relative to SNe II is = 0.12 ( |30|) (Fig. 4). This 

is consistent with observation-based estimates that the 
SNe la frequency is as low as 10 % of total supernova 
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occurrence ©• 

With this relative frequency, ®®Fe from SNe la is about 
50 % of total ^®Fe from all SNe. For the Ch models, the 
abundance ratios between neutron-rich species and ^®Fe 
provide an important constraints on the progenitor sys¬ 
tem. The central density of the white dwarf at thermonu¬ 
clear runaway must be as low as ;S 2 x 10® g cm“^to 
avoid overabundances of ®®Ni, ^"‘Cr, and ^®Ti relative to 
^®Fe, although the exact density depends on the flame 
speed (|30|). Such a low central density can be realized by 
the accretion as fast as M > 1 x IO^^Mq year“^, which 
is consistent with the SD-Ch scenario. SNe la yields thus 
contribute Fe enrichment to the chemical evolution of the 
galaxies. Because the progenitors of SNe la are low-mass 
stars with relatively long lifetimes, Fe could be used as 
a clock to look back into the timing of galaxy formation 

Cosmological Parameters 

Huhble constant: SNe la are very good, but not perfect 
standard candles. The accuracy of the determination 
of Hq using SNe la depends on (i) whether the peak 
absolute-magnitude dispersion of SNe la is sufficiently 
small, and (ii) whether a precise absolute-magnitude cali¬ 
bration of a SN la can be made. Recent progress has been 
made in the empirical determination of Hq. For the dis¬ 
persion, the correlation between the maximum brightness 
and the decline rate (or the light curve shape) is taken 
into account. Multicolor light-curve shapes are also used 
for further corrections (|^. This reduces the dispersion 
from 0.4 to 0.2 magnitudes in the Hubble diagram for 
SNe la (H, 11^) (Fig.5). For the calibration, the distances 
to several host galaxies of SNe la (SNe 1895B, 1937C, 
1960F, 1972E, 1981B, 1989B, and 1990N) have been de¬ 
termined from Hubble Space Telescope observations of 
Cepheids in these galaxies (^). On the basis of these 
data, Hq is estimated to be 58 tg (H); 63.1 ± 6.3 (|3[), 
64 ± 6 (||), and 55 ± 3 ( |^ ) km s“^Mpc“^. 

Successful theoretical models can, in principle, give the 
absolute maximum luminosities, thus providing good es¬ 
timates of Hq and Hm. The maximum luminosity Lmax 
of SNe la when explained with radioactive decay models 
has been used to estimate Hq (|^). From a comparison 
with individual SNe la, Hq = 67 ± 9 km s“^Mpc“^ (p4|), 
and 59 ± 13 km s“^Mpc“^ (13) ■ addition, a fit of the 
late-time synthetic spectra to observations gave Hq = 68 
± 13 km s“^Mpc“^ (1^). 

These estimated values tend to converge to Hq 58 
to 65 km s“^Mpc“^ where the error bars overlap. The 
higher Hq is in conflict with the age of the universe for 
Hm = 1 and Oa= 0, whereas the lower Hq could avoid 
this problem. Further efforts to reduce the error bars will 
require the reexamination of sampling, extinction correc¬ 


tions, and so on. 

Density parameter. In determining and Ha from 
SNe la, absolute calibration of the distances to SNe la is 
not needed, but more accurate measurements of magni¬ 
tudes and z values are necessary. In addition to the dis¬ 
persion problem, whether a significant value of Hm can 
be obtained depends on (i) whether a suitable sample 
of remote SNe la can be obtained, and (ii) whether the 
peak luminosities of SNe la are sufficiently free from the 
effects of cosmic and galactic chemical evolutions. 

More than 60 SNe la with high z values have been 
observed (|3^, |4^) (Fig. 5). Values derived from the first 
7 supernovae at z > 0.35 are Hm = = 0 

or Hm = 1—Ha= 0.94+Q;2g for Hm + Ha =1, assuming 
that evolutionary effects are small. To clarify whether the 
nature of high-z SNe la is the same as that of nearby SNe 
la, it is important to identify the progenitors’ evolution 
and population ( ^ . Systematic studies on the effects of 
metallicity of the progenitors are underway. 
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Figure 1: Spectra of SNe la (SN 1994D and SN 1990N) 
about 1 week before maximum brightness. SN1991T was 
peculiar (|^. The y-axis is magnitude where the units of 
the flux /jy are ergs s“^ cm“^ Hz“^. 

Figure 2: Empirical family of the visual light curves of 
SNe la, which shows the brightness-decline rate relation. 
The triangles, squares, circles, and diamonds denote SNe 
1991T, 1981B, 1986G, and 1991bg, respectively (||). 
Figure 3: The non-LTE spectra of carbon deflagration 
model W7 (unmixed) at 20 days and 23 days past explo¬ 
sion, which are compared with SNe 1992A and 1994D, 
respectively (^^. F\ is the flux in units of ergs s“^ 
Angstrom”^. 

Figure 4: Solar abundance pattern based on synthesized 
heavy elements from a composite of SNe la and SNe H 
with the most probable ratio (30). The dashed lines in¬ 
dicate typical uncertainties of a factor of 2 involved in 
the observational and theoretical abundances, where the 
abundance ratios are normalized to unity for ^®Fe. The 
open and filled circles are used to clarify the association 
to elements. 

Figure 5: Hubble diagram for the seven high-z SNe 
la out to z ~ 0.6, with some of the low-z SNe la (|3|): 
(A) uncorrected blue magnitudes tob, and (B) ms with 
corrections of the brightness-decline rate relation 
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The solid curves are theoretical tob values for (Om, 

= (0, 0) (top curve), (1,0) (middle curve) , and (2,0) 
(bottom curve). The square points are not used in the 
analysis, because these points are corrected based on the 
extrapolation outside the range of light curve widths of 
low-redshift supernovae. 
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